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L.  M.  Woods,1  T.  L.  Reinecke,2  and  A.  K.  Rajagopal3 
1  Department  of  Physics,  University  of  South  Florida,  Tampa,  Florida  33647,  USA 
2Naval  Research  Laboratory,  Washington,  DC  20375,  USA 
3 Center  for  Quantum  Studies,  George  Mason  University,  Fairfax,  Virginia  22030,  USA 
and  Inspire  Institute,  Inc.,  McLean,  Virginia  22101,  USA 
(Received  19  October  2007;  published  29  February  2008) 

The  dynamics  of  spins  in  semiconductor  quantum  dots  often  is  controlled  by  their  hyperfine  coupling  to 
nuclear  spins.  We  develop  a  straightforward  and  efficient  approach  to  describe  the  dynamics  and  the  effective 
decoherence  of  the  electron  spins  due  to  hyperfine  coupling  in  realistic  quantum  dots.  Systems  with  a  large 
number  of  nuclei  and  an  arbitrary  initial  nuclear  polarization  for  which  the  number  of  nuclei  initially  flipped 
over  is  much  less  than  the  total  number  of  nuclei  are  treated.  This  treatment  employs  a  pole  approximation 
within  a  Schrodinger  equation  of  motion  for  the  state  of  the  coupled  electron  and  nuclear  spin  system,  and  it 
allows  us  to  treat  systems  with  arbitrary  initial  conditions.  We  find  that  typical  time  scales  for  the  effective  spin 
decoherence  are  on  the  order  of  tens  of  microseconds. 
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The  spin  degree  of  freedom  of  electrons  is  of  considerable 
interest  in  contexts  ranging  from  “spin  electronic”  devices  to 
implementations  for  quantum  bits  (qubits)  in  quantum  com¬ 
puting.  A  major  obstacle  in  these  applications  is  decoherence 
due  to  coupling  of  the  spin  to  its  environment.  For  example, 
a  spin  in  a  quantum  dot  (QD)  is  an  attractive  candidate  for  a 
qubit,  but  the  spin  decoherence  time  must  be  longer  than  the 
times  needed  to  manipulate  the  spins  in  gate  operations. 

Sources  of  spin  decoherence  in  quantum  dots  include  pho¬ 
non  induced  modulations  of  spin-orbit  coupling1,2  and  hyper¬ 
fine  coupling  to  nuclei.3'4  Recent  work  on  spin  relaxation  due 
to  phonon  scattering  gives  relatively  long  spin  relaxation 
times  at  low  temperatures.1  2  For  sufficiently  high  magnetic 
fields,  electron  spin  lifetimes  on  the  order  of  milliseconds 
have  been  measured.5  At  low  magnetic  fields,  however,  the 
hyperfine  interaction  between  the  electron  spin  and  the  un¬ 
paired  nuclear  spins  typically  present  in  III-V  semiconductor 
quantum  dots  provides  an  unavoidable  source  of  losses  and 
dominates  the  spin  decoherence  in  dots.  The  effective  deco¬ 
herence  in  these  systems  arises  from  variations  of  hyperfine 
coupling  across  the  dot.  Although  strictly  speaking  the 
coupled  electron-nuclear  spin  system  evolves  coherently  in 
time,  coupling  to  a  large  number  of  nuclei  gives  an  effective 
decoherence  in  realistic  laboratory  times.  QDs  have  large 
numbers  of  nuclei  ( 1 04 — 1 06)  with  varying  initial  nuclear  po¬ 
larizations.  Each  of  these  two  features  presents  formidable 
difficulties  in  calculating  the  spin  dynamics  of  these  systems. 

An  approach  for  the  electron  spin  dynamics  for  the  case 
of  a  fully  polarized  nuclear  system  has  been  given,3  and  nu¬ 
merical  simulations  for  models  of  QDs  with  very  small  num¬ 
bers  of  nuclei  {N^  20)  have  been  made. 4,6,7  In  the  case 
where  the  external  magnetic  field  is  much  larger  than  the 
hyperfine  coupling,  the  relaxation  times  can  be  described  by 
a  standard  perturbative  approach.3  A  master  equation  ap¬ 
proach  also  has  been  developed  in  which  the  hyperfine  inter¬ 
action  is  switched  on  adiabatically.9  This  method  is  appli¬ 
cable  to  starting  conditions  for  which  the  electron  and  spin 
systems  are  decoupled,  which  is  generally  not  the  case  for  an 
electron  spin  interacting  with  nuclei  in  solids.  Recently,  an 


analytical  model  based  on  the  Heisenberg  equation  of  motion 
for  high  magnetic  fields  and  for  the  nuclear  interaction 
switched  on  adiabatically  has  also  been  given.10 

Here,  we  present  a  straightforward  method  to  describe  the 
spin  dynamics  in  realistic  QDs  with  large  numbers  of  nuclei 
and  with  varying  initial  polarizations.  It  is  based  on  a 
Schrodinger  equation  of  motion  approach  for  the  wave  func¬ 
tion  describing  the  system  of  coupled  electron  and  nuclear 
spins,  and  it  uses  a  pole  approximation  that  is  applicable  for 
systems  with  large  numbers  of  nuclei.  It  provides  an  efficient 
method  for  calculations,  and  it  allows  us  to  treat  systems 
with  initially  coupled  electronic  and  nuclear  systems. 

We  consider  the  spin  s  of  an  electron  in  the  conduction 
band  in  a  QD  coupled  to  the  nuclear  spins  I,  via  the  hyper¬ 
fine  Fermi  contact  interaction.  For  simplicity  we  take  the 
nuclear  spins  to  have  magnitude  of  1/2.  An  applied  magnetic 
field  B  along  the  z  direction  gives  a  Zeeman  splitting,  and  the 
total  Hamiltonian  is 

H=(e+  y,)iz+  y_s++  y+s_,  (1) 

7*  =  E  AJiv  y_(+)  =  E  zr /,,-(+)  ■  (2) 

i  i  ^ 

Here,  e=  ptBgB  where  ptB  is  the  Bohr  magneton  and  g  is  the 
Zeeman  g  factor.  A,  is  the  hyperfine  coupling  of  the  electron 
spin  in  the  QD  to  a  nucleus  where  A,=Ai>0 ;//( r;) | 2  and  A  is 
the  total  strength  of  the  hyperfine  coupling  in  a  unit  cell. 
u0=  I  / n0  is  the  inverse  density  of  nuclear  spins,  and  iff r,)  is 
the  electron  wave  function  at  the  nuclear  position  r,.  Thus,  A, 
varies  over  the  dot,  and  for  a  dot  with  N  nuclear  spins,  A, 
—AIN.  Dipole-dipole  coupling  becomes  important  only  at 
times  longer  than  10-4  s  and  is  not  included  here.8  The 
nuclear  Zeeman  interaction  is  much  smaller  than  that  for  the 
electron  and  is  not  included. 

The  time  evolution  of  the  electron  spin  is  given  by  the 
state  vector  \<f(t))  of  the  coupled  electron-nuclear  system 
from  ift\ift(t))  =  H\i/4t)),  where  h  is  taken  to  be  1.  The  z  com¬ 
ponent  of  the  total  spin  of  the  coupled  system,  M.=s, 
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+^=]liv  is  a  conserved  quantity,  and  the  total  number  of 
eigenstates  is  2,v+ 1 .  For  a  given  direction  of  the  electron  spin 
and  for  m  of  the  N  nuclear  spins  flipped  down,  there  are 


Nm  = 


arrangements  of  the  nuclear  spins  in  the  QD.  The  total  num¬ 
ber  of  states  for  a  given  M.  is  NM  =Nm+Nm+ 1,  where  we  see 
that  states  with  both  m  and  m  + 1  nuclear  spins  flipped  down 
contribute  to  the  states  of  a  single  Mz.  Details  of  the  classi¬ 
fication  and  enumeration  of  the  states  of  this  system  as  well 
as  their  use  in  finding  the  formal  exact  solution  of  the  time- 
dependent  Schrodinger  equation  are  given  in  Ref.  11.  An 
arbitrary  state  vector  of  the  system  can  be  written  as  a  linear 
superposition  of  orthonormal  states  with  different  Mz, 


I<Am_(,i.)W>  =  2  +  2  X{m+i}(i)|t ;{m  +  1}). 

Nm  Nm+ 1 

(3) 


Here,  |{w})  =  |{JJ-}m{TT},v_m),  with  ra=0, 1 , ...  ,N,  denotes  all 
arrangements  of  the  m  nuclear  spins  that  are  flipped  down 
among  the  total  N  nuclear  spins,  and  (},  J.)  are  the  electron 
spin  states.  Y{m+1}(r)  and  T{m}(?)  are  the  respective  time- 
dependent  probability  amplitudes. 

The  solutions  are  obtained  for  each  value  of  M„  and 
the  complete  solution  is  the  sum  over  all  states  |i fi{t)) 
The  state  is  normalized  giving 


i  -  2  |cv.(m)|‘ 

MJm) 


2  l^{m}(0|  +  2  |^{m+l}(l)|2 


■  (4) 


The  time  dependences  of  the  average  electron  spin  compo¬ 
nents  are 


^2  K 


Mini)  I 


2  \X{m+1](t)\2-^\Y{m](t)\2 

Nm+ 1  Nm 


(5) 


(s+{t)) 


2  2  C M  (m-\)C m , 

M,(m)  Mz(m- 1)  '  Nm 

(6) 


{s_(,t))  =  {sM\  (7) 

where  s+=sx±  isy. 

It  is  straightforward  to  show  that  for  a  given  Mz,  the 
Schrodinger  equation  with  the  Hamiltonian  in  Eq.  (1)  gives 
the  following  closed  set  of  simultaneous  equations  for  the 
amplitudes: 

N  A. 

iY[m}{t)  =  -B{m}Y{m}{t)+  2  ix{m+i }J(t),  (8) 

j*{m}  z 


{m}  A, 

—  7?{m+t}j2f{m+t}j(i)  +  Y  1  |/;‘(t)  •  (9) 


Here,  R{m}=(e+A-S|"’1}A/)/2,  B{m+1]J=B{m}-Aj/2,  and  A 
=  2^jA;.  The  sum  in  Eq.  (8)  is  over  all  nuclei  exclud¬ 

ing  the  ones  flipped  down  for  the  arrangement  {in},  and  the 
sum  in  Eq.  (9)  is  over  all  nuclei  flipped  down  for  the 
arrangement  {to}.  We  have  introduced  the  following  nota¬ 
tions:  X{m+1}j  is  the  amplitude  associated  with  the  electron 
spin  up  and  m  + 1  nuclear  spins  down,  which  is  obtained  by 
flipping  one  of  the  N-m  up  nuclear  spins,  labeled  by  j,  thus 
adding  one  more  nuclear  spin  to  the  previous  //(-down  spins; 
Y{m+\\n  is  the  amplitude  of  a  state  with  to  nuclear  spins  down 
obtained  by  flipping  the  /th  nuclear  down  spin  to  up  in  the 
group  of  to+  1  down  nuclear  spins.  Equations  (8)  and  (9)  are 
converted  to  linear  simultaneous  algebraic  equations  with 
general  initial  conditions  at  /  =  ()  by  taking  Laplace  trans¬ 
forms.  The  exact  solutions  for  the  amplitudes  are  found  as 
follows: 


: 


1 


n  V  (A;/ 2)2 

i(o  +  B,m\-  2j  . 


j*  {»>} 


iw  -  B 


X 


iY[m}{  0)+  2 


iX{m+lll(0)A,/2 


iw  -  B 


{m+1},/ 


{m}  N 

■2  2 

i=l  /=£{/«} 


AiAiY{m+lyi(co)/4 


iw  -  B 


{ra+1},/ 


(10) 


X{m+  1}(W)  = 


1 

iw-B{m+i} 


{in} 

iX{m+ 1}(0)  +  2 

i=l 


y{,n+l}/,(Q)A,/2 
iw  +  B{m+i}/i 


M 


+  2 


A./2 


AlX{m+  !},//,■(  w) 


iw  +  B 


{m+\}/i  l=t{m} 


(ii) 


We  now  consider  the  solutions  of  Eqs.  (10)  and  (11)  in 
time  using  an  inverse  Laplace  transform.  The  large  number 
of  nuclei  precludes  obtaining  exact  results.  However,  we 
note  that  in  the  case  of  predominant  initial  nuclear  polariza¬ 
tion  m<6N  one  can  make  useful  approximations  for  the  am¬ 
plitudes.  The  set  of  equations  for  the  amplitudes  Xrm\(w)  and 
T{m}(o>),  and  subsequently  for  X{m}(r)  and  T{m}(r),  can  actu¬ 
ally  be  given  as  an  infinite  series  of  terms.12  There  are  di¬ 
mensionless  contributions  in  X{m}(f)  and  Yim\(t)  proportional 
to  AjAj  with  denominators  that  are  second  order  in  quantities 
such  as  B{my  For  large  N  and  m  <  N.  these  terms  are  small  to 
the  order  of  1/A2  and  we  drop  them.  In  addition,  there  are 
“self-energy”  corrections  to  quantities  such  as  of  order 
A2  in  the  denominators  of  Xy^w)  and  Yim\(w),  and  we  also 
drop  them.  Here,  we  use  the  fact  that  the  individual  A, 
~  AIN  are  small  compared  to  A,  and  e+A  dominates  when  it 
appears  in  combination  with  any  subset  of  A,  when  m<N.  In 
effect,  we  have  replaced  a  series  of  poles  for  the  amplitudes 
X{m}(w)  and  T{m}(«)  by  a  single  pole,  hence,  the  term  “pole 
approximation”  for  this  approach.  Then,  the  solutions  are 
found  to  be 
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Y{m](t)  «  F{m}(  oyv 

{m}  N 

y'=l  /^{m}  ^)(^{m}  —  ^{m+l}/l) 


AA 


X(1  — B{m+\ }/i)t/2\ 


-X{m+l](0)eiBW‘12 


- — - (1  -e-'(BW+S{-n+l  }.')f/2),  (12) 

/#{m}  4(e  +  A) 


{w+l}  N 


X 


Ai/ 


j=l  l¥={m+l}  e  "*■  '4)(^{m+l},/  B{m+ 1 } j) 


X(1  —  g_'(B{m+l},/_B{m+l}j),/2) 


+  yw(0)e-aWl« 


{/n+1} 


A; 


X  y  — — — -  ei(B{m}+B{m+l}j)ta\  (13) 

j=i  4 (e  +  A) 


We  note  that  the  results  in  Eqs.  (12)  and  (13)  give  the  nor¬ 
malization  in  Eq.  (4)  correctly  for  all  frequencies  to  order  A2 
and  that  these  results  give  the  exact  result  for  the  short  time 
limit  t— >0.  This  approach  gives  accurate  results  for  large  N 
and  m<6N  providing  that  e+A  in  the  denominator  does  not 
become  small. 

For  the  purpose  of  illustrating  these  results,  we  choose  an 
initial  state  in  which  m  spins  are  flipped  down  and  in  which 
each  arrangement  of  the  m  nuclear  spins  is  equally  likely. 
The  initial  probablity  is  taken  to  be  different  for  the  two 
directions  of  the  electron  spin.  This  can  be  thought  of  as  the 
equal  probablity  choice  of  the  starting  state  and  is  specified 
by  the  initial  conditions  Ysm\(t=Qi)=ymHNm  and  Xsm\(t= 0) 
=xj  \INm,  where  ym  and  xm  are  constants  and  all  other  am¬ 
plitudes  are  zero.  We  find  the  following  analytical  expres¬ 
sions  for  the  expectation  values  of  the  electron  spin  compo¬ 
nents: 


(sM  ■ 


~\lK 


!  -  \ym\2] 


{m} 


1  ^4. : 

1-Re—  22r-^T7G{m} 
Nm{m}  j=i  2(e  +  A) 


(14) 


(15) 

^ m  {m} 

Here,  2{,„}  means  summing  over  all  possible  arrangements  of 
the  initially  down  spins.  The  following  notation  is  also  used: 

G{m](t)=  2  -  AP~ - [l-ArW®]. 

/=£{m}  **{m}  “{m+l}/l 

Here,  h  is  recovered  to  give  the  correct  units.  The  function 
Gsm\  is  present  if  there  is  at  least  one  nucleus  flipped  down, 
and  G{,„}  is  zero  if  the  nuclei  are  initially  fully  polarized.  For 
m<N,  the  sum  in  G{m}  is  2)"^,  , ~2jv,  and  it  is  converted 


t  (|ls) 


FIG.  1.  (a)  The  time  dependence  of  (sz)  from  Eq.  (13)  for  initial 
nuclear  polarization  of  75%  is  shown,  (b)  The  time  dependence  of 
(s+)  from  Eq.  (14)  for  initial  nuclear  polarization  of  75%  is  shown. 
We  take  A  =  90  /ueV  (GaAs)  and  e=A. 

into  an  integral  2^— >^/c/3r.  A  similar  treatment  can  be 
made  for  the  summation  in  Eq.  (15).  Inserting  the  definition 
for  and  performing  the  summation,  one  obtains  that 
J-l{m}e2iBf m^~e2i(e+A),lhJ{t)m,  where  J{t)  =  AZ\eiA‘m-  For 
these  evaluations,  the  summations  are  again  turned  into  inte¬ 
grations  here. 

In  Fig.  1,  we  show  the  time  dependences  of  the  electron 
spin  components  for  75%  initial  nuclear  polarization  for  N 
=  104.  For  illustration,  we  have  chosen  a  spherical  quantum 
dot  of  radius  R,  and  use  |  ^(r,)  [2  =  __,n^R/2y  e~4r •  R  for  the  elec¬ 
tron  wave  function.  We  note  that  there  are  rapid  oscillations 
in  ( s+(t )).  They  arise  from  the  exponential  factor  e'(e+A)r/ft 
The  slower  decays  arise  from  the  contributions  of  Aj~AIN 
in  both  ( sz(t ))  and  ( s+(t )).  The  oscillations  in  the  (sz(t))  dy¬ 
namics  are  more  damped  as  compared  to  the  ones  in  (s+(t)) 
due  to  the  prefactor  A//2[B{m}-B{m+1}//].  The  time  depen¬ 
dences  of  ( sz{t ))  and  ( s+(t ))  are  always  nonexponential. 
Strictly  speaking,  the  coupled  system  of  electron  and  nuclear 
spins  is  secular  and  does  not  have  dissipation.  Revivals  of 
the  electron  spin  will  occur  at  much  longer  times  than  those 
shown  here.  Nonexponential  time  decay  was  also  found  in 
other  works.3,9  The  evolution  of  the  decay  envelope  also  has 
been  found  to  be  nonexponential  using  analytical  methods 
for  systems  with  predominant  initial  nuclear  polarization  in 
which  the  nuclear  interaction  is  switched  on  adiabatically10 
and  in  studies  in  which  numerical  approaches  were  used.6,7 

A  limiting  case  of  these  results  is  given  by  that  for  ini- 
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0.0  100.0  200.0 


t  (|is) 


FIG.  2.  The  time  dependence  of  (sz)  from  Eq.  (16)  for  the  fully 
polarized  case  with  hyperfine  coupling  A  =  90  /xeV  and  e=A. 


dally  fully  polarized  nuclei.  For  this  case,  we  take  the  prob¬ 
ability  of  the  initial  electron  being  up  to  be  5c0  and  the  prob¬ 
ability  of  its  being  down  to  be  y0.  Then,  we  find  that 


(sM 


NMtoM.VoI2 


N 

XRe  2 

;=i 


4  (e  +  A)2 


(1  _  e~2i{e+A)tlfi ^iA-tl2h\ 


(16) 


rapid  oscillations  coming  from  contributions  ~  e~2,(€+A^tlfl  and 
slower  envelop  from  contributions  ~e,Aitl2h.  We  note  that  in 
this  case,  the  time-dependent  contributions  are  of  order  of 
~1  IN  smaller  than  (sz(t= 0)),  a  point  that  has  been  made 
earlier.3’9,10  In  the  fully  polarized  case,  the  long  time  enve¬ 
lope  of  (sz(t))  has  a  time  scale  of  80  /xs,  which  is  longer  than 
the  time  scale  for  the  75%  nuclear  polarization  above. 

The  long  time  behavior  of  (sz(t))  and  (,y±(f))  also  can  be 
obtained  by  evaluating  Eqs.  (14)  and  (16)  at  small  distances 
from  the  center  of  the  dot.  In  this  way,  we  find  that  for  the 
initially  partially  polarized  system  with  m  spins  down, 
(sz(t))~  (ft/Af)36,!+0/2  for  long  times  and  that  for  the  initially 
fully  polarized  system  ( sz{t))~(fil At )3/2,  which  is  consistent 
with  the  results  of  Refs.  3  and  9.  We  find  that  the  long  time 
behavior  of  (sz  ±(t))  depends  on  the  dimensionality  of  the 
system  and  on  the  shape  of  the  quantum  dot.3,9  The  explicit 
form  of  the  decay  functions  from  Eqs.  (14)— (16)  allow  one  to 
obtain  results  for  varying  dot  sizes  and  forms. 

In  summary,  we  have  given  a  method  to  calculate  the  time 
dependence  of  electron  spin  components  coupled  by  hyper¬ 
fine  interactions  to  a  large  number  of  nuclei  with  predomi¬ 
nant  initial  polarization  in  quantum  dots.  This  method  is 
straightforward  and  accurate,  it  allows  us  to  treat  efficiently 
systems  with  general  starting  conditions,  and  it  can  be  used 
to  study  relatively  wide  range  of  systems  with  varying  di¬ 
mensionality,  quantum  dot  shape  and  initial  nuclear  spin  po¬ 
larization. 


{sM~\wi{e+A),lh.  (17) 

In  Fig.  2,  the  time  evolution  for  this  case  is  given.  There  are 
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